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ABSTRACT 
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Flapping-wing propulsion is studied experimentally through thrust measurements 
and flow visualization. The objective of the research is to provide further insight into the 
aerodynamics of flapping- wing micro air vehicles (MAVs). Experimental work is 
conducted in the NPS 1.5 m x 1.5 m in-draft wind tunnel. A previously constructed 
model is suspended by thin wires and is used to measure the thrust performance of the 
flapping-wing MAV . For this experiment, the model is tested in four configurations; 
three with varying wing mount stiffness and the fourth with an articulated pitch 
mechanism. Thrust is indirectly determined using a laser range-finder to measure stream- 
wise displacement of the model. Three methods of flow visualization are attempted to 
gain further insight into the flow-field around the MAV. First tufts are placed on and 
around the model to identify the flow-field. Second, a smoke rake placed outside the 
tunnel is used to route smoke into the test section. Thirdly, a smoke wire system is used 
to produce smoke in the test section. Experimental results are compared with flow 
visualization results and previous experimental and numerical work. 
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I. 



INTRODUCTION 



A. OVERVIEW 

In this paper the effects of wing mount stiffness on the aerodynamic performance 
of a flapping-wing Micro Air Vehicle (MAV) was investigated. In addition, flow 
visibility methods were utilized in an attempt to gain further knowledge of the flow field 
about the MAV. While there are many potential propulsive methods for MAVs, 
including conventional or rotary wing designs, this paper builds upon previous research 
conducted at NPS to increase the knowledge base of the flapping-wing propulsion system 
for the MAV application. 

The experimentation for this paper was conducted in the NPS 1.5 m x 1.5 m in- 
draft wind tunnel using a previously developed parallel flapping-wing MAV. To obtain 
thrust data, the MAV was suspended with wire from the roof of the tunnel with thrust 
indirectly measured using a laser rangefinder to measure the streamwise displacement of 
the MAV. 

B. BACKGROUND 

With a shift towards more diverse military operations in recent years, an 
operational need for locally owned and operated reconnaissance assets has been 
identified. The recent development of new technologies such as micro-electromechanical 
systems (MEMS), tiny CCD-array cameras, equally small infrared sensors and chip-sized 
hazardous substance detectors has made the development of small airborne 
reconnaissance vehicles viable. To this end, the Defense Advanced Research Projects 
Agency (DARPA) is funding research into the development of MAVs. A MAV is 
defined by DARPA as “an affordable, fully functional, militarily capable, flight vehicle, 
limited to 15 cm in length, width or height.”[Ref. 1]. These vehicles must be able to stay 
aloft for 20 to 60 minutes, cover a distance of up to 10 km and carry a payload weighing 
up to 20 grams. The flight regime in which MAVs will fly, is for the most part 
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unstudied. The low Reynolds number flight, which will be encountered by MAVs, has 
been up until this point, strictly the domain of insects and small birds. 

Essential to the success of a MAV design is the requirement for high-density 
propulsion and power sources to fulfill the MAV’s propulsive requirements. Previous 
research indicates that a parallel flapping-wing may produce higher propulsive efficiency 
than a conventional propeller driven or rotary winged vehicle in low Reynolds number 
flow. The research summarized in this paper builds upon previous research conducted in 
the study of flapping-wing MAVs. 



C. FLAPPING- WING PROPULSION 

In 1909 and 1912, Knoller [Ref. 2] and Betz [Ref. 3] respectively, independently 
observed that a flapping-wing creates an effective angle of attack, ote, resulting in a 
normal force vector with both lift and thrust components. The Knoller-Betz effect is 
illustrated in Figure 1. 




Figure 1 . Effective Angle of Attack Produced by Flapping-Wing. 
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In 1922, Katzmayr [Ref. 4] demonstrated the Knoller-Betz effect experimentally 
by measuring an average thrust with a stationary wing in a sinusoidally oscillating wind 
stream. This field of research was further advanced in 1924 when Bimbaum [Ref. 5 & 6] 
observed the conditions that lead to flutter or the generation of thrust and suggested a 
flapping-wing as an alternative to the conventional propeller. Over the next few decades, 
much research was conducted on the aerodynamics of flutter and gust response. 

However, it wasn’t until the mid-1930s that further research was conducted by von 
Karman and Burgers [Ref. 7], to explain the theory behind the generation of thrust or 
drag in a flapping-wing. For a more detailed history of flapping-wing research refer to 
Refs. 8 & 10. 

More recently, numerical and experimental work was conducted by Jones and 
Platzer [Ref. 8]. In their research, they found that the two-airfoil, opposed-plunge 
configuration was the configuration which appeared most promising for the MAV 
application due to this configuration’s relatively high propulsive efficiency and the 
inherent stability resulting from the opposed plunge motion. In this application, the 
lower airfoil acts as an image airfoil in a ground-effect analysis. The opposed plunge 



configuration is shown in Figure 2. 
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When operating in ground effect, downwash velocities are decreased. The 
ground presents a boundary condition requiring the normal component of velocity to go 
to zero altering the streamline pattern around the airfoil. In an analytic ground-effect 
analysis, a mirror image of the airfoil is placed across the plane of symmetry. The vortex 
system of this image airfoil is opposite the original system. In order to satisfy the 
boundary condition, the vertical velocities induced by these two image vortex systems 
will cancel at the plane of symmetry. In addition to vertical velocity going to zero, the 
destructive interference of the counter-rotating vortex systems diminishes the influence of 
the trailing vortices on the airfoil. 

The parallel flapping-wing configuration studied in this paper approximates the 
beneficial aspects of ground-effect when flying out of ground-effect. In this three- 
dimensional case, the reduced influence of the tip vortices, lessen the induced angle of 
attack, increasing lift and thrust. 
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II. EXPERIMENTAL ANALYSIS 



A. EXPERIMENTAL APPARATUS 
1. MAV Construction 

The MAVs used in this research were constructed by Dr Kevin Jones [Ref. 10]. 

A side view of the MAV with the articulated pitch mechanism can be seen in Figure 3. 
The MAV is designed around an exceptionally small, geared stepping motor from RMB 
Smoovy. The motors are 5mm in diameter with the 25:1 planetary gear system, are 
approximately 25mm long and weigh a mere 2.4g. These motors produce about 2.5 xlO' 3 
Nm of torque at speeds up to 800 RPM. 




Figure 3. MAV with Articulated Pitch. 



The brushless motors are controlled by an oscillatory driver circuit, which allows 
for very precise, steady and reproducible rotational speeds. As the control circuit and 
power supply are large and bulky, they are currently mounted external to the model. 
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Three-pole power is fed to the model through four 0.076mm diameter copper wires which 
support the model. The support wires connect to the model via small gold plated pins, at 
the nose and near the tail of the model. The mass of the wires is negligible compared to 
the mass of the model and they are flexible enough so as not to impede the model’s 
motion. 

The model’s body has composite construction, built primarily out of balsa wood 
and very thin graphite-epoxy laminates. The motor is connected by thin-wall silicon 
tubing to a crankshaft that moves the flapping beams via Scotch yokes constructed of 
piano wire. The wing’s leading edge is a teardrop shaped piece of balsa wood with 
carbon-fiber ribs running to the trailing edge of the Japanese tissue wing surface. In this 
paper, the wings had a span of 150mm and a chord of 35mm. 

The wings are super-glued to the flapping beams using thin carbon fiber strips 
whose length and width control the elasticity of the joint. The elasticity of the wing 
mount allows for a passive feathering mechanism. The wing deflects in pitch 
proportionally to the moment about its leading edge. In this paper, the elasticity of the 
wing joint was adjusted by decreasing the distance between the wing’s leading edge and 
the flapping beam. Three different distances were used for the tested model. In addition, 
an active pitch control was added to the model by gluing arms and strings to the model’s 
body and wings to achieve a pure plunge motion for the wing. The four configurations 
were tested in the NPS low-speed in-draft smoke tunnel and the results compared. 

2. Wind Tunnel 

Experiments were performed in the NPS low-speed in-draft smoke tunnel 
depicted in Figure 4. The tunnel was modeled after the smoke tunnel of the Naval Air 
Engineering Lab in Philadelphia [Ref. 12]. Air is ingested from inside the building 
through a square 4.5 x 4.5 m inlet, converging through a 9: 1 bell shaped contraction to a 
1.5 x 1.5 m test section. Tunnel speed is controlled by a variable pitch fan driven by a 
constant speed electric motor. Motor and fan vibration are isolated from the test section 
by rubber sleeves on each side of the motor/fan assembly. The tunnel velocity range is 0 
to 9.5 m/s. 
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-igure 4. NPS Low-Speed In-draft Wind Tunnel. 



Late in 1999, Lt Tim Lund rebuilt the neglected wind tunnel [Ref. 1 1]. His 
extensive repairs included fixing the fan, resetting the variable pitch fan blades to achieve 
lower velocities, removed old fixtures, refinished the tunnel floor and the old pitot-static 
system was replaced. A tunnel calibration was performed using LDV following tunnel 
repairs and was used in the reduction of the data for this paper. In addition, the test 
section turbulence intensity level was determined and remained below 1.75% with tunnel 
velocities over 1.5 m/s [Ref. 11]. 

3. Miscellaneous Equipment 

A pitot-static tube located approximately two feet below and one foot upstream of 
the model was used to determine the tunnel velocity. The pitot-static tube was connected 
to a MKS Baratron type 223B differential pressure transducer. The transducer voltage 
output, which is linear with pressure, was read by a Metex 3610D digital multimeter that 
linked the voltage data to a PC for recording. 
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A frequency strobe light was used to set the flapping-wing frequency. The strobe 
was set to the desired frequency, and the wing frequency adjusted to match. The 
accuracy of the strobe was checked using a Monarch Instrument Tach IV Digital Optical 
Tachometer. Strobe error varied with frequency to a maximum of 1.5 percent in the 
frequencies of interest. Actual frequency used for data reduction was determined after 
the experiment using voltage data from the model controller. 

A Matsushita Electric Works Model ANL1651AC laser rangefinder was used to 
determine model displacement when flapping. The sensor accuracy as prescribed by the 
manufacturer is ± 100pm ± 0.002 x Ax for the range 130mm ± 35mm. Its measurable 
range is 8 - 18 cm, with accuracy decreasing with increased distance from the center 
point. The rangefinder was positioned behind the model at its center point distance of 13 
cm. During runs, this corresponded to a position 3 to 4 chord lengths behind the model 
and was not thought to create a significant flow interference effect. The output of the 
rangefinder is 1 V/cm displacement. 

A Link Instruments DSO-2102 digital storage oscilloscope was used to measure 
voltage. The DSO-2102 was connected to a desktop PC for display and recording of the 
voltage signals. The oscilloscope simultaneously records two data channels. One 
channel was designated for the laser rangefinder, and this data was subsequently used to 
determine thrust. The second channel was used to collect controller voltage data for 
eventual determination of flapping frequency and to synchronize data. The oscilloscope 
features 32 kilobyte buffers per channel. 

A Setra EL410D digital balance was used to measure the mass of the model 
before each run. The scale has an accuracy of ± O.OOlg. The model’s mass was found to 
vary by as much as 0.005g from day to day, due to humidity, dust or other environmental 
contamination. A Rosco model 4500 smoke generator was used for flow visibility. A 
smoke rake was constructed to match the output of the smoke generator. 

A Sony DCR-VX1000 digital Video camera was used to capture video images of 
the MAV in motion. Using the manual focus, auto exposure, a shutter speed of 1/8 sec 
and the strobe light set at approximately 12 Hz (just slightly faster than the MAV 
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flapping frequency), acceptable video images of the MAV were obtained. Capturing 
flow visualization images using the smoke rake or smoke wire proved to be more 
difficult. A Kodak DC 260 digital camera was used to capture flow visualization images 
(thread tell-tales). Best results were achieved by setting the manual focus to 18 inches, 
the exposure to 0.5 sec (the slowest manual timed exposure) and the strobe light and 
MAV flapping frequency set to approximately 8 Hz. This configuration lead to some 
blurring as approximately 4 cycles were captured per image. 

B. EXPERIMENTAL CONFIGURATIONS 



Wing Span 


b* 


150 mm 


Wing Chord 


C 


35 mm 


Aspect Ratio 


AR=b 2 /S=b/c 


4.29 


Taper Ratio 


A,=c t /c 0 


1 


Sweep Angle 


A 


0° 


Velocity 


v„ 


0-6 m/s 


Flapping Frequency 


F 


8, 10, 12, 14 Hz (approximate) 


Mean Wing Separation 


Ay m 


0.8429c 


Plunge Amplitude 


H 


25.4 mm 


Reduced Frequency 


K 


0.29 - °o 



Table 1. General Experimental Configuration. 



Wing Stiffness (Low) 


K 


4.746 x 10" Nm 


Wing Stiffness (Medium) 


K 


6.384 x 10' 3 Nm 


Wing Stiffness (High) 


K 


14.17 x lO' 3 Nm 


Wing Stiffness (Articulated Pitch) 


K 


N/A 



Table 2. Wing Mount Stiffness. 
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c. 



EXPERIMENTAL METHOD 



Before each wind tunnel run, the model with mounting pins was weighed using 
the Setra scale. The model was then hung in the wind tunnel with four 0.076mm 
diameter copper wires and mounting pins as can be seen in Figure 5. After the model 
was installed in the tunnel the laser rangefinder was positioned near its zero voltage 
center point, 13 cm behind the model. The DSO program was then started on the PC and 
setup to record data at 1 KHz (approximately 32 sec of data). The model was given a 
push, and once the displacement fell within the limits (500 mV/division, 10 divisions) 
displayed in the DSO program (± 2.5 cm), data recording began. Three separate sets of 
data were taken before either the tunnel or the model was turned on. This data was used 
to determine both the effective pendulum length for thrust calculations and the model’s 
initial position. A close-up image of the MAV in motion with laser can be seen in Figure 
6 . 




Figure 5. MAV Installed in Tunnel. 
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Figure 6. MAY in Motion and Laser. 



Before any thrust measurements were taken, a 60-second sample of pitot-static 
data (sample rate of 1 Hz) was taken to establish a zero velocity voltage. This data was 
recorded on the PC using the software provided with the Metex 3610D DMM. 

Next, a set of thrust measurements were taken with the tunnel off. With the 
controller mode selector knob in the variable speed position, the power was turned on. 
The default power setting for this start-up mode was set so that the model’s wings would 
be flapping at approximately 8 Hz. The strobe was set to 8 Hz to verify the flapping 
frequency. In order to record frequency data, the TTL output from the controller was 
connected to the DSO. With the 25:1 gearbox installed in the model, 25 voltage pulses 
from the controller correspond to one complete cycle for the model. Once the model 
position stabilized, recording of the position data by the DSO program, set to sample at 
10 KHz, was started. Once again, three sets of data were recorded at each frequency. 

The strobe frequency was then increased in 2 Hz increments up to 14 Hz. At each 
2 Hz increment, the signal to the MAV was stepped up until the flapping frequency 
matched the strobe or the motor stalled. The motor tended to stall if the load became too 
great or the MAV was disturbed by an external disturbance. 
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Once the static measurements were complete, the wind tunnel was turned on. 
Initially, the pitch on the fan blades was set to give zero tunnel velocity. The pitch was 
slowly increased incrementally to the velocity of interest. Once the velocity stabilized, 
60-seconds of pitot-static data was recorded on the PC (1 Hz sample rate). 

As the focus of this paper is on determining the feasibility of using a flapping- 
wing mechanism to generate thrust, the drag of the MAV is required so that the gross 
thrust generated by the MAV can be determined. As the drag of the MAV varies with 
wing position (maximum drag when wings are together), the wings were flapped at an 
extremely low frequency, on the order of 0. 1 Hz, to determine an average drag. The DSO 
was used to record the position data from the laser rangefinder at 1 KHz sample rate. 
Again three sets of data were taken. The position data from the DSO for the stiff wing 
mount run at approximately 5 m/s and 0.1 Hz can be seen in Figure 7. 



Stiff Wing Mount Model - Raw Drag Data 
Flapping Frequency - 0.1 1 Hz 
Tunnel Velocity - 5.14 m/s 




Sample 

Figure 7. Change in MAV Drag with Wing Position. 

Once the drag data had been recorded, thrust measurements were taken as per the 
static case at 8, 10, 12 and 14 Hz (if the motor allowed). Once the measurements were 
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made, the tunnel speed was increased until a tunnel speed of 5 to 6 m/s was reached. At 
tunnel speeds in this range, the model’s drag displaced it sufficiently to force the 
displacement data outside the display limits of the DSO software. With tunnel speeds 
from 0 - 5.9 m/s and flapping frequencies from 8 - 14 Hz reduced frequencies in the 
range k = 0.29 - °° were observed. 

D. EXPERIMENTAL DATA REDUCTION 

At each test point, three different sets of data were recorded; voltage from the 
differential pressure transducer, output from the controller and voltage data from the laser 
rangefinder. To reduce the pressure transducer data to velocity [m/s]. Equation 1, 
obtained from an LDV- calibration completed by Lt Lund [Ref. 1 1], was used. 

Vel = 0.4819(raV) 0 489 (!) 

The pressure data file was opened in Excel and converted to velocities using 
Equation 1. The average and standard deviation of the 60-second sample were then 
taken. To account for changing environmental conditions, an initial and final pressure 
transducer reading were taken. A linear change in the zero condition was assumed and 
the difference between the initial and final condition was divided by the number of 
velocities and added to each velocity reading. The test velocities ranged from 0 - 5.9 
m/s. To determine the model’s effective pendulum length for thrust calculations and the 
model’s initial position, the 32-second (1 Hz sample rate) static swing data file was 
opened in Excel and the model’s displacement was plotted versus sample number (see 
Figure 8). 
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Figure 8. MAV Displacement vs Sample #. 

Next, the data was edited so that an integer number of cycles were represented. 
The number of cycles and the number of samples within the cropped data set were 
counted and the frequency of the hanging model was determined. The installed model’s 
frequency,^, was approximately 1.07 Hz for the configurations studied. The model’s 
pendulum length, L, was then calculated using Equation 2 where g is the acceleration due 
to gravity. 




The model’s initial position was determined by averaging the cropped data set. 
Next a MATLAB program (Appendix A), based upon one written by Lt Lund [Ref. 11], 
was written to reduce the raw voltage data from the rangefinder and the controller and to 
complete error calculations. To convert the voltage data to model displacement in 
centimeters, d , Equation 3 was used: 



Medium Stiff Wing Mount Model 
Static Swing Data 
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d = vpd *5.12 — 



Vj*vpd* 10.24 
256 



( 3 ) 



Where vpd is the volts per division setting in the DSO software display, V is the 
voltage output from the laser rangefinder, d\ is the initial displacement of the model in 
centimeters, 10.24 is the number of divisions displayed and 256 is the number of pixels 
displayed. 

The force required to displace the model was calculated using Equation 4: 



Where M is the mass of the model. Equation 4 was used to determine both the MAV’s 



force required to displace the MAV model was subtracted from the net thrust. 

As both the thrust and drag of the MAV vary throughout its flapping cycle, the 
data reduction program was required to ensure an integer number of cycles were used for 
its calculations. To accomplish this, the voltage data from the controller was analyzed 
and cropped as it passed through a certain trigger voltage. One voltage cycle is 
equivalent to a motor revolution from one positive peak to the next positive peak. One 
can see from Figure 9 that the controller gives a crisp, repeatable voltage spike for each 
cycle. As the gearbox has a 25:1 reduction ratio, 25 peaks equate to 1 flapping cycle for 
the MAV. After the MATLAB program identifies the start of the cycle, it counts the 
number of times the controller voltage passes through the trigger voltage (on the positive 
slope) until the 25 th time when the program calculates the voltage, thrust, and frequency 
for the flapping cycle. 

As three samples of data were taken at each velocity/frequency setting, the means 
from each sample set were averaged. The reduced data was saved to a summary file for 
compilation upon completion of the program. The program also allows for a 
comprehensive summary file for each velocity/frequency setting. This comprehensive 
summary was useful for trouble shooting during construction of the program. 



F = Mg * tan arcsin 

v 1 )j 




(4) 



drag and its net thrust at a certain tunnel speed. In order to obtain gross thrust, the drag 
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Controller Voltage [Pixels] 



Controller Voltage 




Sample # 

Figure 9. Controller Voltage. 
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E. EXPERIMENTAL RESULTS AND ANALYSIS 



1. Experimental Results 

In this paper, the effects of wing mount stiffness on the performance of the 
flapping-wing MAV were investigated. Four configurations were tested; three with 
increasing wing mount stiffness and one with an articulated pitch mechanism which 
attempted to obtain pure airfoil plunge. 

In Figures 10 through 13, the thrust obtained from each configuration is plotted 
versus velocity at approximately 8, 10, 12 and 14 Hz flapping frequency. For three- 
dimensional plots of the same data, see Appendix B. 



Low Stiffness Model 




Figure 10. Thrust Performance of Low Stiffness Model. 
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Thrust [N] 



Medium Stiffness Model 




Figure 11. Thrust Performance of Medium Stiffness Model. 

High Stiffness Model 




Figure 12. Thrust Performance High Stiffness Model. 
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Articulated Pitch 




Figure 13. Thrust Performance of Articulated Pitch Model. 

From Figures 10 through 13 one can see that the general shape of the thrust vs 
velocity plot at different frequencies is basically the same for each configuration. 

Initially, the model produces a high thrust, which quickly reduces to a minimum as 
velocity increases. When the tunnel speed is further increased, the thrust developed by 
the model begins to slowly recover. One can see that as the flapping frequency increases, 
the thrust- velocity lines move upward and to the right for each configuration. 

Figures 14 through 17 compare the thrust performance of the four different 
configurations at each frequency. One can see that in the static case and at extremely low 
tunnel velocities, the models with the more flexible wing mounts produce greater thrust. 
The articulated pitch model does not perform well in the static and low speed flow cases. 
One can see that as the tunnel velocity increases up to the velocity which produces 
minimum thrust, the thrust performance of the low and medium stiffness models and the 
high stiffness and articulated pitch models approach each other. Once past the velocity 
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Thrust IN] 



that produces minimum thrust the performance of the medium stiffness and articulated 
pitch models surpass that of the high stiffness and low stiffness models. 



Frequency - 8 Hz 




Figure 14. Configuration Thrust Comparison @ 8 Hz. 
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Frequency - 10 Hz 




|~-+— Low Stiff — Medium Stiff — x— High Stiff — x— Articulated PitcTT 

Figure 15. Configuration Thrust Comparison @ 10 Hz. 



Frequency - 12 Hz 




Figure 16. Configuration Thrust Comparison @ 12 Hz. 
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Frequency - 14 Hz 




Figure 17. Configuration Thrust Comparison @ 14 Hz. 

From the above results, it appears that the best wing mount configuration for a 
MAV of similar construction to the one used in these tests would be one that had a 
flexible mount at low velocities and a high stiffness or articulated pitch (pure plunge) 
mount at higher speeds. 

2. Error Analysis 

There are four primary sources of error associated with the experimental thrust 
measurements. These four include the error due to the average deviation of thrust 
between cycles, laser rangefinder error, oscilloscope error, and weight/wire length error. 
The sum of these four errors is presented as the vertical error bars in Figures 10 through 
13. As the in-draft tunnel vents to the atmosphere, the velocity measurements taken by 
the pitot-static system were highly sensitive to wind conditions. Figures 10 through 13 
also illustrate the fact that the velocity error varies greatly day to day. On the days that 
the data for the low stiffness model was tested, the wind was relatively calm. However, 
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the second day that the medium stiffness model was tested was extremely windy (> 
30mph). The resulting error in velocity was as much as 22%. The velocity error bars 
displayed in Figures 10 through 13 represent one standard deviation added and subtracted 
from the mean tunnel velocity. 

The first source of thrust error is thrust average deviation. The mean thrust for 
every complete flapping cycle of the model was calculated. This cycle mean thrust of the 
model varied from cycle to cycle. This variation was due to a variety of factors but the 
most significant contributors were that the tunnel intake was indoors and that the outlet 



vents to the atmosphere. Although the tunnel room was vented to the outside by several 
banks of windows, the opening and closing of doors and other windows in the building 
caused slight velocity drift which affected model displacement and thrust. In addition, 
variations in wind velocity also affected results. Average deviation is similar to standard 
deviation in that it provides a measure of the width of the data about the mean. It is a 
more robust estimator, and is given by Equation 5 [Ref. 1]. 



1 N 

ADev(x 1 ...x N ) = — £ K " *| 



(5) 



j = i 



For this application, x represents the mean of each cycle, j is the cycle counter, N is the 
total number of cycles, and x is the mean of all cycles. The average deviation was 
calculated at every frequency and velocity and is considered here to be an estimated error 
bound. It is the predominant source of error for velocities above approximately 4 m/s. 

The second source of thrust error was the laser rangefinder. Rangefinder error 
consists of a fixed error at the center point (13 cm). A linear error proportional to the 
distance from the center point is added to the fixed error at all other points. For the 
ranges encountered here, the rangefinder accuracy is given by the Equation 6: 

Range Error = 0. 1 mm ± 0.2% * Distance from Center Point (6) 

The third source of thrust error was due to the digital oscilloscope. An error due to 
the maximum resolution of the oscilloscope was added to all voltage measurements. This 
error varies with the selected maximum value of the displayed vertical axis. The highest 
accuracy is obtained when the voltage signal approaches the maximum displayable, but 
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does not exceed the display window. The oscilloscope accuracy analysis indicated that 
the actual oscilloscope error was always less than half of its resolution. Statistically, this 
is reasonable in that errors would tend to occur on both sides of the actual voltage, 
offsetting each other. In the interest of a conservative analysis, however, the maximum 
resolution was used as the error. Due to the requirement to capture the maximum 
displacement of the model at tunnel velocities approaching 6 m/s, the DSO resolution had 
to be set at 0.5 Volts per division. Consequently, the error associated with DSO 
resolution was the greatest contributor to thrust error. 

The fourth source of error was that associated with the thrust calculations. As 
discussed in the Experimental Data Reduction section, to determine the thrust produced 
by the model the weight and length of the model’s wires had to be determined. In order 
to account for the uncertainties in these measurements, a thrust error calculation was 
carried out using the maximum weight and minimum wire length and was subtracted 
from the mean thrust for each data set. 
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III. FLOW VISUALIZATION 



A. FLOW VISUALIZATION METHODS 

Three methods of flow visibility were used to try and gain insight into the flow 
field about the MAV. Initially, a smoke rake and fog machine located outside the tunnel 
intake were used to feed smoke into the test section. Then tell-tale tufts were placed on 
and around the MAV’s wings. Finally a smoke wire smoke system was used to bum fog 
fluid in the test section near the MAV. 

1. Smoke Rake and Fog Machine 

Initially, attempts to document the flow field around the MAV were conducted 
with a smoke rake located just outside the tunnel intake connected to a Rosco 4500 
smoke generator by a 25 foot long six-inch ID flexible duct. A smoke rake was 
constructed with an eight foot-long main tube of six-inch ID, schedule 40 PVC pipe with 
16-1.5 inch ID smoke tubes cemented to the main tube. This rake was mounted on a 
height adjustable stand and can be seen in Figure 18. 




Figure 18. Smoke Rake In Front of Tunnel Inlet. 
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After several attempts a couple of problems with using the smoke rake located in 
front of the tunnel intake were discovered. As the fog fluid is burned to generate smoke, 
it is much warmer than ambient air. Consequently, when the smoke enters the bottom of 
the rake, it tends to jet past the first few smoke tubes and most of the smoke exits the rake 
from the top eight tubes. The smoke continued to rise as it entered the tunnel’s 9: 1 
contraction. 

In addition, as the smoke passed through the contraction, the smoke stream tubes 
began to wrap around each other, forming a large spiral of smoke by the time it reached 
the test section. As the wing span of the MAV was only 150mm and the diameter of the 
smoke spiral was approximately 30 cm in diameter at the test section, the smoke totally 
enveloped the MAV and this method was fairly useless for flow visibility purposes. 
Consequently, other methods of flow visibility were investigated. 

2. Tell-Tale Tufts 

As Jones [Ref. 10] had some success using single strands of cotton thread (a 
thread is usually constructed of three or more intertwined strands) glued to the MAV’s 
wings, this method was attempted. To reduce the strand inertia and stiffness problems 
associated with tufts glued to the wings, single strands of cotton were glued to a balsa 
wood spar constructed from a MAV wing’s leading edge. Even though the strands were 
not attached to the rapidly flapping-wing, problems with the inertia and stiffness of the 
fibers were encountered. Once an optimal fiber length was determined (approximately 
40 mm) the spar was placed in various positions around the MAV. The best results were 
obtained with the spar glued to the center rear of the MAV’s fuselage with the strands 
passing between the two flapping-wings. 

3. Smoke Wire System 

The third method used to document the flow field about the MAV was a smoke 
wire system. For this method, a bare strand of 0.022” nickel chromium (nichrome) wire 
was fixed to span the tunnel just ahead of the MAV in the test section. Electrical 
connectors were attached to the end of the wire on the floor of the tunnel, and on an 
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adjustable brass tube just above the MAV. The adjustable brass tube also had a port to 
which a clear silicone tube was attached which delivered fog fluid from a pressurized jug 
outside the test section. The jug was pressurized using a bike pump through one of two 
adapters bolted through the jug’s cap. The clear silicone tube was attached to the second 
adapter and fog fluid flow was controlled using a pair of hemostat clamps. 

Power for the smoke wire was provided by a Variac set to provide 15-20V. The 
Variac setting was dependent upon tunnel speed and fog fluid flow rate. The model was 
installed in the tunnel on a sting that was mounted on a 60 cm span symmetric airfoil that 
was in tum mounted to a two axis traverse. The flow was illuminated using a Smith- 
Victor Corp 750W Photographic Light. 

Different configurations of the wire were tried to determine the optimum design. 
Initially, loops were bent in the wire so that the fog fluid would pool there and produce 
smoke for an extended period. This method did cause the fluid to pool but the loops 
created significant vortex shedding. Next, 0.005” bare copper wire was tied around the 
nichrome wire to create longer duration smoke. Vortex shedding was once again 
observed, however the magnitude was much less than that caused by the looped nichrome 
wire. The looped wire worked fairly well, however after a few fog cycles, the copper 
wire knots slid down the nichrome wire. Finally, a straight bare nichrome wire was used. 
Though some vortex shedding was still observed at higher flow speeds, this configuration 
gave the most consistent smoke plane. The smoke wire set-up can be seen in Figures 19 
and 20. 
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Figure 19. Smoke Wire Set-Up. 



Figure 2U. ^moke Wire bet-Up m tunnel. 
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Figure 21. Tell-Tale Tufts, Medium-Stiff Wing, 12Hz Static. 

With the given frequency and exposure time, approximately six cycles were 
captured in this image. This can clearly be seen in the multiple images of the fourth 



Flow images from the smoke wire technique were captured using a Sony DVC- 
2000 digital video camera and a Kodak DC 260 digital camera. When using the Kodak 
camera, the focus was manually set at either 24 or 36 inches and the flash was turned on 



B. FLOW VISIBILITY RESULTS 



1. Tell -Tale Tufts 

Using the tuft method, the best results were obtained when the tuft spar was 
placed just ahead of the leading-edge plane of the wings mounted on the MAV fuselage. 
In Figure 21, the medium-stiff winged MAV was flapping at 12 Hz with no wind. The 
MAV was lit by a strobe light synchronized with the model and the image captured using 
the Kodak DC 260 digital camera with the focus set manually to 18 inches and the 
exposure set to 0.5 sec. 
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strand from the right. These multiple images of this strand indicate that the flow between 
the wings, near the centerline of the MAV is quite unsteady. This image shows the three 
dimensionality of the flow field around the MAV as the flow is entrained as the wings 
move apart. Figure 22 illustrates the outward flow from between the wings as they move 
together. This image is of the MAV with the medium-stiff wing mount, flapping at 12 



Hz with flow speed of approximately 2.9 m/s. 




Figure 22. Tell-Tale Tufts, Medium-Stiff Wing, 12Hz, 2.9 m/s Wind Speed. 



2. Smoke Wire System 

Initially, attempts were made to illuminate the smoke flow using the strobe light. 

It quickly became clear that the output of the strobe would not be sufficient for capturing 
smoke images. In addition, as there were no means available to synchronize the camera 
with the strobe, an exposure length longer than one complete cycle was necessary to 
ensure that the strobe flash was captured. As the flow around the MAV is unsteady, these 
attempts to capture an image over more than one cycle lead to significant blurring of the 
smoke plane. Finally, a Smith-Victor Photographic Light was used to illuminate the flow 
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and the Kodak DC260’s flash was turned on. The images in Figures 23 through 26 were 
taken from approximately three to four feet downstream of the MAV with the manual 
focus set to two or three feet. 




Figure 23. MAV Wing Tip Vortex Induced Wake Structure. 



In Figure 23, the smoke wire is positioned approximately six inches upstream of 
the MAV’s port leading edge. One can clearly see the wake structure developed by the 
wing tip vortices as they are shed during the flapping motion. Looking upstream, as the 
upper and lower wings move apart, they respectively shed counter-clockwise and 
clockwise rotating vortices. Correspondingly, when the upper and lower wings move 
together, they respectively shed clockwise and counter-clockwise rotating vortices. The 
wing tip vortices can be more clearly observed in Figure 24. 
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Figure 24. Close Up of MAY Wing Tip Vortices. 



The vortices that are clearly visible in Figure 24, were formed as the upper and 
lower wings moved together. As these vortices grew and coalesced, they formed a jet, 
which flowed outwards away from the MAV’s centerline. The interaction of these 
vortices also pushed the tip vortices outwards. Once the wings reached their minimum 
separation, they stopped shedding their vortices and the smoke plane remained flat as can 
be seen in the vertical smoke line just outboard of the wing tips joining the two vortices. 
When this image was taken, the MAV’s wings had just begun to move apart. The tip 
vortices which were just forming, are not visible but the effects of the inward flowing jet 
they formed, can be seen in the deflected smoke plane between the wings to the right of 
the smoke wire in the image. 

The alternating vortices of the upper wing and their effects on the smoke plane 
can be seen in Figure 25 taken from above and downstream of the MAV. In this image 
the clockwise rotating vortex nearest the MAV was shed from the upper wing as the 
wings moved together. The effects of the outward flowing jet created by this vortex set 
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on the upper portion of the smoke plain can be seen in the deflection of the top of the 
plain to the left (seen to the left of the vortex in this image). 



Figure 25. Alternating Vortices Shed from the MAV’s Upper Wing. 

The second vortex aft of the MAV is counter-clockwise rotating and was shed 
from the upper wing as the wings moved apart. The effects of the inward flowing jet 
created by this second vortex set on the upper portion of the smoke plain can be seen in 
the deflection of the top of the plain to the right (seen above and to the left of this vortex 
in the image). 
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The final image obtained by the smoke wire method can be seen in Figure 26. 

This image was taken with the smoke wire positioned in front of the MAV’s nose directly 
along its centerline. In this image one can see the effects of the boundary layer over the 
MAV’s nose in the upward deflection of the smoke plane. As this image was taken, the 
wings were at their minimum separation; consequently, the motion of the wings is 
approaching zero and the smoke plane is relatively undisturbed above the upper wing. 
Just behind the upper wing, the smoke plane is displaced upwards by the clockwise 
rotating (as seen from this side of the smoke plane) vortex shed from the upper wing as 
the wings moved together. As expected in Figure 26, the flow is relatively two 
dimensional when compared to the flow about the MAV’s wing tips. 




Figure 26. Two Dimensional Flow Over the MAV’s Centerline. 



Figures 27 through 35 summarize of the flow field information obtained using the 
smoke wire system. In Figure 27, the MAV’s wings are at their minimum separation and 
the vortices shed during the closing of the wings have been convected downstream. 
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Figure 27. MAV Flow Field - Wings at Min. Separation - Steady. 

In Figure 28, the MAV’s wings are beginning to move apart. The wings are 
deflected due to aerodynamic loading and the wings’ inertia. Due to the movement of the 
wings and the resulting change in lift generated by this movement, vortices are formed 
off the wings’ trailing edges and wing tips. 
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"igure 28. MAV Flow Field - Wings at Min. Separation - Moving. 
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In Figure 29 the wings have moved apart and accelerated to their maximum 
plunge speed at their median separation. As the wing tip vortices grow they begin to 
coalesce and form a jet between them, flowing towards the MAV’s centerline. The 
interaction of the vortices also pushes the vortex pair towards the MAV’s centerline. The 
wind convects the vortices downstream. 
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Wind Direction 



Figure 29. MAV Flow Field - Wings at Median Separation. 

In Figure 30 the wings are decelerating as they reach their maximum separation. 
The intensity of the shedding votices decreases and they continue to be convected 
downstream. 
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'igure 30. MAV Flow Field - Wings Approaching Max. Separation. 

In Figure 31 the MAVs wings are at their maximum separation and are no longer 
moving. Consequently, the wings cease to shed their vortices. 
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'igure 31. MAV Flow Field - Wings at Max. Separation - Steady. 
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In Figure 32 the MAV’s wings are beginning to move together. One can see that 
the wings deflect in the opposite direction and the vortices are rotating in the opposite 
sense to those in Figure 28. 
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figure 32. MAV Flow Field - Wings at Max. Separation - Moving. 



In Figure 33 the wings have moved apart and accelerated to their maximum 



plunge speed at their median separation. As the wing tip vortices grow they begin to 
coalesce and form a jet between them, flowing away from the MAV’s centerline. The 



interaction of the vortices also pushes the vortex pair away from the MAV’s centerline. 
The wind convects the vortices downstream. 
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In Figure 34 the wings are decelerating as they reach their minimum separation. 
The intensity of the shedding votices decreases and they continue to be convected 
downstream. 
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*igure 34. MAV Flow Field - Wings Approaching Min. Separation. 



39 




Finally, in Figure 35 the MAVs wings are at their minimum separation and are no 
longer moving. Consequently, the wings cease to shed their vortices, which have been 
convected downstream. 

MAV in 
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Wind Direction 

Figure 35. MAV Flow Field - Wings at Min. Separation - Steady. 
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IV. CONCLUSIONS 



With a shift towards more diverse military operations in recent years, an 
operational need for locally owned and operated reconnaissance assets has been 
identified. The recent development of new technologies has made the construction of 
small airborne reconnaissance vehicles viable. Nature hints that flapping wings are 
efficient at the low Reynolds numbers in which these MAVs would operate. 
Consequently, the potential efficiency of flapping wings as a means of MAV propulsion 
makes investigation worthwhile. 

Experimental testing of a previously constructed MAV was conducted in the NPS 
1.5 m x 1.5 m in-draft wind tunnel. For this paper, four different wing configurations 
were investigated; three with varying degrees of wing mount stiffness and one with an 
articulated pitch mechanism which attempted to obtain a pure-plunge motion for the 
wings. Thrust was indirectly determined by measuring the stream-wise displacement of 
the suspended model due to flapping. The steady-state drag was subtracted from the 
computed thrust in order to isolate the thrust performance of the MAV and minimize 
viscous boundary-layer drag effects. 

After collection, reduction and analysis of experimental data, it became apparent 
that in order to optimize the thrust performance of a MAV some sort of active pitch 
control would be beneficial. At low flow speeds, a flexible wing would produce more 
thrust, while at higher speeds a stiffer mount or active pitch control to approach a pure 
plunge motion would be more effective. 

Images collected using the smoke wire method were extremely beneficial in 
determining characteristics of the flow field about the MAV. Further flow visibility 
efforts need to be made at different flapping frequencies and tunnel velocities so that the 
flow mechanics responsible for the changing thrust performance can be better 
understood. 
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V. 



RECOMMENDATIONS 



In order to move forward in the development of the flapping-wing as a propulsion 
system for MAVs, more effort must be put into visualization of the flow field about the 
MAV. Initial attempts using the smoke wire system were fairly effective at capturing 
images, but the method can be improved. First, a constant temperature circuit employed 
in the smoke wire system would simplify the smoke generating process. Secondly, the 
use of a laser sheet would allow images of the cross section of the wake behind the MAV 
to be captured. Thirdly, a higher intensity strobe light or multiple synchronized strobes 
may allow photographic equipment to capture “still” images of the flowfield. This would 
allow the photographer more control over which portion of the cycle is captured on film. 
Finally, a flow highlighting system originating from the MAV wing would allow for 
further investigation into the dynamic stall phenomenon occurring on the MAV wings 
throughout the flapping process and also highlight the formation of the wing tip vortices. 
Perhaps the use of titanium tetra-chloride in the in-draft tunnel could be investigated. 
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APPENDIX A. MATLAB PROGRAM FOR DSO-2102 OSCILLOSCOPE DATA 

REDUCTION 



%Capt Sean Duggan KAY / DSO Data Reduction Program 
%Data Analysis Process: 

%- Weigh model with pins before each set of runs (enter weight Sc max weight [N] below) 

%- Install model in tunnel, center laser, swing model to get period data (1 KHz sample 
rate) 

%- Reduce period data to get wire length [cm], enter below (L= (g) / (2*pi*freq) ~2) 

%- Use period data to determine initial model displacement [pixels] , enter below 
% 

% Record initial velocity (pressure) data before tunnel on 
% Record velocity at beginning of each run 
%- Record final velocity with tunnel off 
%NOTE - True velocity=0 . 4319* (millivolts) '0 . 489 

% (From LDV / Pitot system calibration from Lt Tim Lund's Thesis) 

% 

% IMPORTANT : Save data in following format: 

% - freqi IOC) , vel (vOOO) , data setd thru 3) 

% (ie. lOGvOGGl, 100vG002, 100v00G3, lOlvOOQl, . . . , lxxvG003 etc. for each vel) 

% use freq 100 for drag data (IKHz sample rate) 

% - use freq 101... lxx for thrust data (IGKHz sample rate) 

% - take 3 sets of data at each freq/ vel setting 

% - save data at each velocity in its own directory (ie. vOOO, vOQl, v002...) 

% 

% Archive data on CD-R, go to computer with MATLAB to reduce data 

% After data is reduced at each velocity, copy data in Mean_Thrust.wri file to summary 
file 
% 

%NOTES: If you open Mean_Thrust.wri in EXCEL, copy data to summary file then close 

% Mean_7hrust.wri. If you don't close the file, this program will crash when 

% it tries to delete Mean_Thrust.wri file 

% 

%*- When data reduction complete, subtract drag (IGOvXXXX) from thrust at each vel to get 
% corrected thrust 

%- Compile velocity data in one summary file 

%- Determine avg vel and standard deviation of vel and copy to thrust summary file 



clear 
clear all 



delete H:\docs\Results\Mean_Thrust.wri; 



% ********************* Fill in this part ****************************** 
f ilelocation= ' F : \ ; % CHANGE FOR EACH COMPUTER CD Rom drive location 

date= ' 17 May 00 V ; % CHANGE FOR EACH RUN Date directory 

dimame= ' v018\ ' ; % CHANGE FOR EACH VELOCITY ■ Velocity Directory 

f ilename= ' vOlS ' ; 



filetype=' .pm' ; 
start filenumber =100 ,- 
lastf ilenumber=104 ; 



% First file number - Drag data (0.1Hz flap) 

% CHANGE FOR EACH VELOCITY - Last file number 



voltsperdiv=0 .5; 
weight=0. 21329; 
wei ght_max= 0.21334; 
wire_length=21 . 8087386; 
wire_length_min=21 . 8090208 ? 
disp_intial_pix=123 . 0946956 ; 



% Model weight in Newtons 
% Max weight for error calcs 
% Cable length in cm 
% Min wire length for error calcs 
% Initial displacement at vO fO in pixels 



disp_initial= (voltsperdiv*5 . 12 ) - (disp_intial_j?ix) * (voltsperdiv*10 .24/256) ; 



*■***■*•*•******* 






*:**** + **■# 



f i dl = f open ( ' H : \doc s \Re suits \Mean_Thrus t . wr i ' , ' w ' ) ? 
fprintf (fidl, ' FilenameX tVel\tVel\tMean\tCorrected\ tTotal ErrorXt' ) ; 
fprintf (fidl, ' Avg Dev\tLaser\tDSO\tWt Wire \tKean\fc Avg Dev\tMin\ tMax\n ' ) ; 
fprintf (fidl, '\t\tStd Dev\tThrust\tThrust\t+/ \ tError\tError\tError\t ' ) ; 
fprintf (fidl, ' Error \t?requency\ tFrecuency\t Frequency \tFrequency\n' ) ; 
fprintf (fidl, ' \t [m/s] \t [m/s] \t [Nj \t [N] \t [N] \t [N] \t [N] \t [N] \t [N] \t ' ) ; 
fprintf (fidl, ' [Hz] \t [HzJ \t [Hz] \t [Hz] \n\n' ) ; 
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fclose { f idl ) ; 



for i=startfilenuinber : lastf ilenumber 
for k=l : 3 

clear data 
if i==100 

s amp 1 e_r a t e = 1 0 0 0 ; 
else sample_rate=10000 ; 
end 



Clears data matrix each iteration 



Sample rate in Hz 
Samole rate in Hz 



fm DSO (0.1Hz flap 
fm DSO (Thrust) 



Drag) 



9; ^***** + **** + ***-*-** + v + * Loads data *********x*******^x*>*******-*** 
disp{ ( ' ******* Computing num2str { i) , filename, num2str (k) , filetype, ' ******•*']) 

rawdata=load ( ( f ilelocation, date, dirname, num2str ( i ) , f ilename,num2str (k) , filetype] ) ; 
if k ==1 

datal { : , 1) =rawdata( : ,2) ; % CH 3 is loaded to datal col 1 
datal ( : , 8 ) =rawdata { : , 1 ) ; % CH A is loaded to datal col 8 
% ******** Converts data to voltage and distance (lV=lcm) ********** 
datal { : , 2 ) = (voltsperdiv*5 . 12) - ( (datal ( : , 1 ) ) * ( vol tsperdiv*10 . 24/256 ) ) ; 
datal ( : , 9 ) = (voltsperdiv*5 . 12) - ( (datal ( : , 8 ) ) * (voltsperdiv*10 . 24/256 ) ) ; 

% ******** Does Nothing *********************x*******+*+*x**** 
datal ( : , 3 ) =datal ( : , 2 ) ; % Does nothing 

% ******** Converts voltage to Force (N] ************************** 
datal ( : , 5)= weight * tan (as in ( (datal ( : , 3) -disp_initial ) /wire_length) ) ; 

% ******** Computes the avg dev of thrust * free in data *********** 

% DSO data variables 
if i==100 

sig_trip = 132; 
else sig_trip = 140; 
end 

sig_oldl = 257; 
sig_old2 = 257; 
sig_old3 = 257; 

% Voltage variables 
cycle_volts=0; 
eye 1 e_me an_vo 1 1 s = 0 ; 
volts_suml=0 ; 

% Thrust variables 

cycle_thrust=0; % Sum of thrust for one cycle 

cycle_mean_thrust=0 ; % Mean thrust of one cycle 
thrust_suml=0; % Sum of mean thrust of all cycles 

% Frequency variables 

% Mean frequency of one cycle 
% Sum of frequencies of all cycles 
% Variable for min frequency 
% Kin free counter 



% 



% 



Kin data value from controller (G.l Hz flap - Drag run) 
Min data value from controller (Thrust runs) 

sig_oldl, 2, and 3 are the 3 data points prior 
to the current pt. Use 257 if the trip signal 
is neg voltage. Use -1 if trip is a + voltage. 



% Sum of voltage for one cycle 
> Mean volts of one cycle 
Sum of voltages of all cycles 



cycle_mean_freq=0; 
f req__suml = 0; 
f reql=0 ; 
x=l ; 

% Counter variables 
nr_cyclesl = -l ; 
cycle_data_pts=0; 
total_data_ptsl = 0 ; 

y=-i; 



% Counter for total number of cycles in data 
% Number of data pts -with in one cycle 
% Number of data pts of all cycles 

% Counter to determine cycles (25 trips = 1 cycle) 
data_pointsl=length ( datal (:, 1 )) ; % Counts rows of data matrix 
for j=l :data_pointsl % Looks at all data points 
if abs (datal ( j , 2) ) > 5 

error {'** Voltage is outside rangefinder limits (+/- 5V) **') 

end 

if nr_cyclesl >= 0 % Ignores partial data before 1st cycle 

cycle_volts = cycle_volts + datal ( j -1 , 2) ; % Sums cycle volts 
cycle_thrust = cycle_thrust + datal (j-1, 5) ; % Sums cycle thrust 
cycle_data_pts=cycle_data_pts+l ; % Counts data pts in cycle 

end 

% Is microswitch tripped? 

if datal ( j , 8) >=sig_trip & sig_oldl<sig_trip & ... 

sig_old2<sig_trip & sig_old3<sig_trip % This assumes a neg 



% trip voltage. Change 

y=y+i; 

if nr_cyclesl == -1 
nr_cyclesl = 0; 
end 



>, < signs if trip is a + voltage. 
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if y == 25 

nr_cyclesl=nr_cyclesl+l; % 25 Pulses/cycle 
y=0; % Resets counter 

% Volts 

cycle_mean_volts=cycle_volts/cycle_data_pts; % Cycle mean 
datal { j-1, 4) =cycle__mean_volts; % Stores cycle rr.ean volts 
volts_suml=volts_suml+cycle_mean_vol ts; % Sums total volts 
% Thrust 

cycle_mean_thrust=cycle_thrust/cycle_data_pts; % Cycle rr.ean 
datal {j-1, 6) =cycle_mean_thrust ; % Stores cycle mean thrust 
thrust_suml=thrust_suml+cycle_mean_thrust ; % Sums tot thrust 
% Frequency 

cycle_mean_f req=sample_rate/cycle_data_pts; % Cycle mean 
datal { j -1 , 7 ) -cycle_mean_f req; % Stores cycle mean frequency 
f req_suml = f req_suml+cycle_mean_f req; % Sums all cycle freqs 
if datal {j-1, 7) >0 

freql{x) = cycle_mean_f req; % Save nonzero freqs 
x=x+l ; 

end 

total_data_ptsl=total_data_ptsl+cycle_data_pts; 

% Reset variables 

cycle_volts=0; % Reset cycle -voltage for new cycle 

cycle_thrust=0; % Reset thrust summation for new cycle 
cycle_freq=0; % Reset free summation for new cycle 

cycle_data_pts=0; % Reset nr of data pts for new cycle 

end 

end 

sig_oldl = datal (j, 8); % Last signal value 

if j >=3 

sig_old2 = datal { j -1 , 8) ; % The data value before the last one 
sig_old3 - datal { j -2 , 8) ; % Three data points back 

end 

end 

end 



if k = -2 

data2 { : , 1) =rawdata( : , 2) ; % CH 3 is loaded to datal col 1 

data2 ( : , 8) =rawdata( : , 1) ; % CH A is loaded to datal col 8 

% ******** Converts data to voltage and distance (lV=lcm) ********** 

data2 { : , 2) = ( { voltsperdiv*5 - 12 ) - < (data2 ( : , 1) ) * (voltsperdiv*10 . 24/256 ) ) ) , 

data2 { : , 9) = ( {voltsperdiv*5 . 12 ) - { (data2 ( : , 8) ) * (voltsperdiv*10 . 24/256 ) ) ) , 

% ** + -*•**•** dqqs No thine ************************************** 

data2 ( : , 3 ) =data2 ( : , 2 ) ; % Does nothing 

% ******** Converts voltage to thrust [N] ************************** 
data2 { : , 5) = weight* tan (asin ( (data2 ( : , 3 ) -disp_initial ) /wire_length) ) ; 

% ******** Computes the avg dev of thrust * f req in data *********** 

% DSO data variables 



if i==100 

sig^trip = 132; 
else sig_trip = 140; 
end 

sig_oldl = 257; 
sig_old2 = 257; 
sig_old3 = 257; 

% Voltage variables 
cycle_volts=0; 
eye 1 e_mean_vo lts=0; 
vol ts_sum2=0 ; 

% Thrust variables 
cycle_thrust=0; 



% Min data value from controller (0.1 Hz flap - Drag run) 
% Kin data value from controller (Thrust runs) 

% sig_oldl, 2, and 3 are the 3 data points prior 
% to the current pt . Use 257 if the trip signal 
% is neg voltage. Use -1 if trip is a + voltage. 

% Sum of voltage for one cycle 
% Mean voles of one cycle 
p 6 Sum of voltages of all cycles 



% Sum of thrust for one cycle 
cycle_mean_thrust=0 ; % Mean thrust of one cycle 
thrus t_sum2 = 0 % Sum of mean thrust of all cycles 
% Frequency variables 

cycle_mean_f req=0 ; % Mean frequency of one cycle 

freq_sum2=0; % Sum of frequencies of all cycles 

freq2=0; % Variable for min frequency 

x=l; % Min free counter 
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o Counter variables 

nr_cycles2 = -l ; 
cycle_data_pts=0 ; 
total_data_pts2=0 ; 

y=-i; 



1 Counter for total nurtber of cycle 
% Number of data pts within ore cycle 
% Number of data pts of all cycles 
% Counter to determine cycles (25 



data_points2=length (data2 ( : , 1 ) ) ; °6 Counts rows data mat 
for j = l : data__points2 % Looks at all data points 
if abs (data2 ( j , 2) ) > 5 

error ('** Voltage is outside rangefinder limits { + / 



end 



ir. data 




5V) **') 



cycle 



if nr_cycles2 >= 0 % Ignores partial data before 1st cycle 

cycle_volts = cycle_volts + data2 ( j -1, 2) ; % Sums cycle volts 
cycle_thrust = cycle_thrust + data2 ( j-1, 5) ; % Sums cycle thrust 
cycle__data_pts=cycle_data_pts+l; % Counts data pts in cycle 

end 

% Is micros'll tcfc tripped? 

if data2 ( j , 8) >=sig_trip & sig__oldl<sig_trip & ... 

sig_old2<sig_trip & sig_old3<sig_trip % This assumes a neg 
% trip voltage. Change >, < signs if tr_p is a + voltage. 

y-y+i; 

if nr_cycles2 == -1 
nr_cycles2 = 0; 

end 



if y == 25 

nr_cycles2=nr_cycles2+l; % 25 Pulses/ cycle 
y=0; % Resets counter 

% Volts 

cycle_mean_volts=cycle_volts/cycle_data_pts; % Cycle mean 
data2 ( j -1 , 4 ) =cycle_mean_volts; % Stores cycle mean volts 
volts_sum2=volts_sum2+cycle_mean_volts; % Sums total volts 
% Thrust 



cycle_mean_thrust=cycle_thrust/cycle_data_pts; % Cycle mean 
data2 ( j -1 , 6 ) =cycle_mean_ thrust ; % Stores cycle mean thrust 
thrust_sum2 = thrust_sum2+cycle_mean_thrust; % Sums tot thrust 
% Frequency 

cycle_mean_freq=sample_rate/cycle_data_pts; % Cycle mean 
data2 ( j -1 , 7 ) =cycle_mean_f req; % Stores cycle mean frequency 
f req_sum2 = f req_sum2+cycle_mean_f req; % Sums all cycle freqs 
if data2 ( j-1, 7) >0 

freq2(x) = cycle_mean_f req; % Save nonzero freqs 
x=x+ 1 ; 

end 

total_data_pts2=total_data_pts2+cycle_data_pts; 

% Reset variables 

cycle_volts=0; % Reset cycle voltage for new cycle 

cycle_thrust=0; % Reset thrust summation for new cycle 
cycle_f req=0; % Reset free summation for new cycle 

cycle_data_pts=0 ; % Reset nr of data pts for new cycle 

end 

end 



sig_oldl = data2(j,8); % Last signal value 

if j>=3 

sig_old2 = data2 ( j-1, 8) ; % The data value before the last one 
sig_old3 = data2 ( j-2, 8) ; % Three data points back 

end 



end 



end 



if k == 3 

data3 ( : , 1 ) =rawdata ( : , 2 ) ; % CH 3 is loaded to datal col 1 
data3 ( : , 8) =rawdata ( : , 1) ; % CH A is loaded to datal col 8 
% ******** Converts data to voltage and distance (lV=lcm) ********** 
data3 ( : , 2 ) = ( ( voltsperdiv* 5 . 12) - ( (data3 ( : , 1 ) ) * ( voltsperdiv*10 . 24/256 ) ) ) ; 
data3 ( : , 9) = ( ( voltsperdiv* 5 . 12) - ( (data3 ( : , 8 ) ) * ( voltsperdiv* 10 . 24/256) ) ) ; 

ft Qq 0 g Jlo thiro ******** *■*■*•** + *’>»*••***■**■*••■»•*'*••*'>«*• ■*••*+••*** 

data3 ( : , 3 ) =data3 ( : , 2) ; % Does nothing 

% ******** Converts voltage to thrust [N] ********-**•*************** 
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data3 (: ( 5)= weight* tan (as in ( (data3 ( : , 3) -disp_ini tial ) /wire_length) ) ; 
% --*'-**' Computes the aver dev of thrust & free _r. data -**-**->*-* 
% DSO data variables 
if i==100 



sig_trip = 132; 
else sig_trip = 140; 
end 

sig_oldl = 257; 
sig_old2 = 257; 
sig_old3 = 257; 

% Voltage variables 
cycle_vol ts=0 ; 
eye 1 e_me an_vo 1 1 s = 0 ; 
volts_sum3 =0 ; 



% Kir. data value from controller (0.1 Hr flap - Drag run) 
% Kin data value from controller (Thrust runs) 

‘o sig_oidl, 2, and 3 are the 3 data points prior 
% to the current pt. Use 257 if the trip signal 
% is neg voltage. Use I if trip _s a + voltage. 

% Sun. of voltage for one cycle 
% Kean volts of one cycle 
% Sum of voltages of all cycles 



% Thrust variables 



cycle_thrust=0 ; % Sum. of thrust for one cycle 

cycle_mean_thrust=0; % Kean thrust of one cycle 
thrust_sum3=0; % Sum of mean thrust of all cycles 

% Frequency variables 

cycle_mean_f req=0 ; % Kean frequency of one cycle 

f req_sum3=0 ; % Sum. of frequencies of all cycles 

freq3=0; % Variable for min frequency 



x=l ; % Kin freq counter 

% Counter variables 



nr_cycles3=-l; 
cycle_data_pts=0 ; 
total_data_pt s3 = 0 ; 
y=-l; 



% Counter for total number of cycles in data 
% Number of data pts within one cycle 
% Number of data pts of all cycles 

Counter to determine cycles (25 trips * 1 cycle) 



data_points3=length(data3 ( : , 1) ) ; % Counts rows of data matrix 
for j=l :data_points3 % Looks at all data points 
if abs (data3 ( j , 2) ) > 5 

error ('** Voltage is outside rangefinder limits (+/- 5V) 



end 

if nr_cycles3 >= 0 % Ignores partial data before 1st cycle 

cycle_volts = cycle_volts + data3 ( j -1 , 2) ; % Suras cycle volts 
cycle_thrust = cycle_thrust + data3 ( j-1 , 5) ; % Sums cycle thrust 
cycle_data_pts=cycle_data_pts+l; % Counts data pts in cycle 



end 



% Is microswitch, tripped? 

if data3 (j , 8) >=sig_trip & sig_oldl<sig_trip & ... 

sig_old2<sig_trip & sig_old3<sig_trip % This assumes a neg 
% trip voltage. Change >, < signs if trip is a + voltage. 

y=y+i; 

if nr_cycles3 == -1 
nr_cycles3 = Oz- 
ena 



if y == 25 

nr_cycles3=nr_cycles3+l; % 25 Pulses /cycle 
y=0; % Resets counter 

% Volts 



cycle_mean_volts=cycle_volts/cycle_data_pts; % Cycle mean 
data3 ( j-l, 4) =cycle_mean_volts ; % Stores cycle mean volts 
volt s_sum3 =vol t s_sum3 +cyc le_mean_vol ts ; % Sums total volts 
% Thrust 



cycle_mean_thrust=cycle_thrust/cycle_data_pts; % Cycle mean 
data3 ( j -1 , 6) =cycle_mean_thrust ; % Stores cycle mean thrust 
thrust_sum3=thrust_sum3 +cycle_mean_thrust ; % Sums tot thrust 
% Frequency 

cycle_mean_f req=sample_rate/cycle_data_pts; % Cycle mean 
data3 ( j-l, 7) = eye le_mean_ freq; % Stores cycle mean frequency 
f req_sum3=f req_sum3+cycle_mean_f req; % Sums all cycle freqs 
if data3 ( j-l,7)>0 

freq3 (x) - cycle_mean_f req; % Save nonzero freqs 
x=x+l ; 
end 

total_data_pts3=total_data_pts3+cycle_data_pts ; 

%> Reset variables 
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cycle_volts=0 ; % Reset cycle voltage tor new cycle 

cycle_thrust=0 ; % Reset thrust summation for new cycle 

cycle_freq=0; % Reset free summation for new cycle 

cycle_data_pts=0 ; % Reset nr of data pts for new cycle 

end 

end 

sig_oldl = data3(j,8); % Last signal value 

if j >=3 

sig_old2 = data3 ( j -1, 8) ; % The data value before the last one 
sig__old3 = data3 ( j -2, 8) ; % Three data points back 

end 

end 

end 

end 

% Voltage and frequency calculations: 
data= [datal;data2;data3 ] ; 

total_data_pts=total_data_ptsl+total_data_pts2+total_data_pts3 ; 

vol ts_sum=vol ts_suml+vol ts_sum2+vol ts_sum3 ; 

f req_sum=f req_suml+f req_sum2+f req__sum3 ; 

thrus t_sum= thrus t_suml + thrus t_sum2 + thrus t_sum3 ; 

nr _cyc 1 e s =nr_c y c 1 e s 1 + nr_cy c 1 e s 2 + nr_cyc 1 e s 3 ; 

freq=[freql freq2 freq3]; 

clear datal 

clear data2 

clear data3 

volts_mean=volts_sum/nr_cycles ; 

if volts_mean==0 % No or very slow flap case, still need data 

vol ts_mean=mean (data ( : , 3) ) ; 

end 

inner_suml=0; % Inner summation for adev of volts 
inner_sum2 = 0; % Inner summation for adev- of free 
freq_mean=freq_sum/nr_cycles; % Kean freq of cycles 
max_freq=max ( f req) ; % Max cycle freq of all cycles 
min_freq=min ( f req) ; % Min cycle freq of all cycles 
data_points=length (data ( : , 1 ) ) ; 
for m=l :data_points 
if data(m,6)~=0 

inner_suml = inner_suml + abs (data (m, 4 ) -volts_mean) ; 
inner_sum2 = inner_sum2 + abs (data (m, 7 ) -f req_mean) ; 

end 

end 

adev_freq=inner_sum2/nr_cycles ; % Avg dev of frequency 
% Thrust 

thrust_mean=weight* tan (asin ( ( volts_mean-disp_initial ) /wire_length) ) ; % Mean thrust o 
cycles 

% *** There are 4 sources of thrust error: 

% *** 1st one is due to cyclic nature of flapping wing thrust: ****** 
adev_volts=inner_suml/nr_cycles; % adev fm voltage data 
% *** 2nd source of error laser rangefinder error [mm] ■** + ** + ** + ** 
if abs (volts_mean) < 3.5 

error_laser = 0.1 + 0.002 * abs(10 * volts_mean) ; 
else 

error_laser = 0.25 + 0.005 * abs(10 * volts_mean) ; 

end 

error_laser = error_laser/10 ; % Convert range error to volts error 

% * * * Oscilloscope volt max error=resolution of DSO- 2102 3rd error * 
error_DS0 = (voltsperdiv*5 . 12 ) /128 ; 

% *** Convert errors fm voltage to thrust [Nj error *+**+**+**+**+** 

adev_thrustl=weight* tan (asin ( (adev_volts) /wire_length) ) ; 
error_thrust2=weight*tan (asin ( (error_laser ) /wire_length) ) ; 
error_thrust3=weight* tan (asin ( (error_DSO) /wire_length) ) ; 

% * * T Computes error due to weight and wire length 4th error: -***•<**-*** 

error_thrust4 =we ight_max* tan (asin ( (volts_mean-disp_initial) /wire_length_min) ) - 
thrust_mean; 

error_thrust=adev__thrustl+error_thrust2+error_ thrust 3 +error_thrust4 ; 
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% «*+***"**"*-**+* Saves ail data to fziensme.wri ******v****^****** 

% Used for trouble shooting the program 

%£ ilesavez- { ' 'H: \docs\Resuits\ ' , numlstr ( i ) , filename, ' .v;ri ' ] } ; 

% t icz = f open ( f i lesavez , ' w ' } ; 

%fprintf { f icz , [numzstr ( i } , filename] ) ; 

%rprintf {£id2, ' \rCH 3\ tCH B\ tCH B\ tCH 3\tCH B\t')? 

%fprir*tf (fi<52, 'CK 3\tCK 3\tCH A\ tCH A\n'}; 

%fprintf (rid2 , 'Data\tVoltage\ tAhsol Volt\tCycle Mean \t Thrust \t ' } ; 

% fprinti ( f id2 , 'Cycle Mean\ tFreq Kean \ cDat a \ t Vo 1 1 age \ r. ' ) ? 

% fprintf ( f ic2 , ' \t fV] , [cm] \t [V: f [cm] \t [V] \t [Nj \t[N] \t ' ) ; 

%f print f ( f id2 , ' [Hz ] \t\ t [V] \n ' } ; 

%fprintf { f id2 , ' \r.\tMean thrust =\t\t\t %10 . Sf ' , thrust_mean > ; 

% f pr in t f ( f id2 , ' \ n\ tThrus t Avg Dev Error = \ t \ t \ t % 1 0 . 8 f ' , adev_ thrust 1 ) ; 

%fprintf ( f id2 , ' \n\tLaser Error = \ t \ t \ t %10 . 3f ' , error_thrus~2 > ? 

%fp rintf ( f id2 , ' \n\tOscilloscope Error =\t\t\t %10 . 3f ' , error_thrust3) ; 

%fprintf (f id2, ' \r.\tWt Wire Error =\t\t\t %10 . Sf ' , error_thrust4 ) ; 

%fprintf ( f id2 , ' \n\tTotal Error = \ t \ t \ t %1D . Sf ' , error_thrust ) ; 

% f print f ( f id2 , ' \rl\ t \ t\ t \ tMean frequency = \ t \ t %4 . 2f \r. ' , . . . 

% f r eq_mean ) ; 

% f pr in t f ( f id2 , ' \ t \ t \ t \ t Avg Dev = \ t \ t % 4 . 2 f \n ' , acev_f req ) ; 

%fprintf ( fic2 , ' \t\t\t\ zKin Freq = \t\t %4 . 2f \n ' , min_f req) ? 

%fprmtf { f id2 , ' \ t\t\ t\CMax Freq =\t\t %4 . 2f \n\n' ,max_frec) ? 

%fprintf { f id2 , ' %3.0f\t %4.3f\t %5.4f\t %5.4f\t %8.5f\t %8.5f\t %5.2f\t %3.0f\t 
%4 . 3f \n' , data. ' ) ? 

%f close ( f id2 ) ; 

% + ***** + ** + **•> s aves summary of data to Mean_Thrust.wri *********** 
fid3 = fopen( 'H: \docs\Results\Mean_Thrust.wri 'a' ) ; 
fprintf ( f id3 , [num2str ( i ) , filename] ) ; 

fprintf (fid3, ' \t\t\z%10 . Sf \ t\t%10 . 3f \z%10 . 8f \ t%10 . 8f \t%10 . Sf \ t%10 . 8f \ t%4 . 2f \ t%4 . 2f \t%4 . 2f 
\t%4 . 2f \n' , . . . 

thrust_mean, error_thrust , adev_thrustl , error_thrust2 , error_thrust3 , . . . 
error_thrust4 , f req_mean, adev_f req,min_f req,max_freq) ; 
f close ( f id3 ) ; 

end 

disp ( [ ' ********* Program Complete ******** ** *'] ) 
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APPENDIX B. THREE DIMENSIONAL PLOTS OF EXPERIMENTAL DATA 



Low Stiff Model x 10 3 




Velocity (m/s) 

Figure 36. 3-D Plot of Low Stiffness Model’s Thrust Performance. 



Low Stiff Model 
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Figure 37. Contour Plot of Low Stiffness Model’s Thrust Performance. 



53 




Medium Stiff Model 
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Figure 38. 3-D Plot of Medium Stiffness Model’s Thrust Performance. 




Velocity (m/s) 

Figure 39. Contour Plot of Medium Stiffness Model’s Thrust Performance. 
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High Stiff Model 



xIO ' 2 




Velocity <m/s) 

Figure 40. 3-D Plot of High Stiffness Model’s Thrust Performance. 



High Stiff Mode! x10 2 
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Figure 41. Contour Plot of Medium Stiffness Model’s Thrust Performance. 
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Articulated Pitch Model 




Frequency {Hz) 



Velocity (m/s) 

Figure 42. 3-D Plot of Articulated Pitch Model’s Thrust Performance. 




Velocity (m/s) 

Figure 43. Contour Plot of Medium Stiffness Model’s Thrust Performance. 
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